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ABSTRACT 

We derive quantitative predictions of the optical and X-ray luminosity functions for quasars in the 
redshift range z ;$ 3. Based on accreting black holes as primary sources for the quasar outputs, we 
investigate how the accretion is controlled by the surrounding structures, as these grow hierarchically 
from the formation of the host galaxies to their assemblage into poor and eventually into rich groups. 
We argue that for z < 3 efficient black hole fueling is triggered by the encounters of a gas-rich host with 
its companions in a group; these destabilize the gas and induce accretion, giving rise to the following 
features. The dispersion of the dynamical parameters in the encounters produces luminosity functions 
with the shape of a double power-law. Strong luminosity evolution is produced as these encounters 
deplete the gas supply in the host; an additional, milder density evolution obtains since the interactions 
become progressively rarer as the groups grow richer but less dense. We carry out these arguments to 
derive a specific model for the evolving luminosity functions. From the agreement with the optical and 
the X-ray data, we conclude that the evolution of the bright quasars is driven by the development of 
cosmic structures in two ways. Earlier than z ~ 3 the gas-rich protogalaxies grow by merging, which also 
induces parallel growth of central holes accreting at their full Eddington rates. In the later era of group 
assemblage the host encounters with companions drive onto already existing holes further but meager 
accretion; these consume the gas reservoirs in the hosts, and cause supply-limited emissions which are 
intermittent, go progressively sub-Eddington and peter out. Then other fueling processes occurring in 
the field come to the foreground; we specifically discuss the faint emissions, especially noticeable in 
X-rays, which are expected when hosts in the field cannibalize satellite galaxies with their meager gas 
contents. 

Subject headings: Quasars: general - galaxies: nuclei - galaxies: interactions - galaxies: clusters: 
general - X-rays: galaxies 



1. INTRODUCTION 

We are motivated by the steep evolution of the bright 
quasar (QS) population; out to z w 2.5 their number ex- 
ceeds the local value by some 10 2 (Schmidt & Green 1983, 
Schmidt 1989). Around z w 3 the evolution culminates, 
and farther out it goes into a decline (Osmer 1982, Shaver 
et al. 1996). 

Our work will be based on two widely shared notions. 
First, the QSs are powered by massive black holes (BHs) 
accreting gas from their host galaxies; thus e.m. outputs 

L ~ T]c 2 Am/At i5 10 48 erg s _1 (1.1) 

obtain when the baryonic mass Am is accreted over the 
time At with conversion efficiency up to r\ ~ 10 _1 (Rees 
1984). Second, the structures surrounding the BHs de- 
velop after the hierarchical clustering (see Peebles 1993); 
this envisages the galaxies as baryonic cores within heavier 
halos of cold dark matter (DM). Before z ~ 2.5 galactic 
halos are built up through merging of smaller units; there- 
after the galaxies begin to assemble into small groups of 
dynamical mass Mq <L 51O 12 M0, which in turn merge 
into richer groups and eventually into clusters. By such 
a hierarchical development, the numbers of structures in 
each mass range undergo a fast increase followed by a slow 
demise. 

On these grounds, the high QS luminosities in the op- 
tical, IR or X-ray bands may be expected to signal out to 
z ~ 5 environmental conditions conducive to accretion. 



In fact, the decline (a rise in cosmic epoch) of the 
early QSs is widely held to be geared to galaxy forma- 
tion. Specifically, their environment is held to be con- 
stituted by protogalactic spheroids being built up hierar- 
chically through major, chaotic merging events between 
dense, very gas- rich subunits. In these events the gravita- 
tional potential is strongly distorted, while the baryons are 
lumped and shocked; so the angular momentum j of the 
gas orbiting at kpc distances in the host is efficiently trans- 
ferred to the heavier DM substructures, and the gas can 
begin an inward course toward pc or smaller scales. From 
the plenty of gas so made available, central massive BHs 
collapse and/or accrete rapidly; they can accrete at their 
full self-limiting rate, and easily attain Eddington lumi- 
nosities Le ~ 10 47 Mbh/W 9 Mq ergs -1 . This is indicated 
by eq. (1.1), when the BHs double their mass and make 
Am ~ Mbh over short dynamical scales At that match 
the effective Eddington time r\tE ~ 51CU 2 Gyr. In sum, 
the number of early QSs ought to track the galaxy forma- 
tion, and rise together with the number of host galactic 
halos exceeding some 1O 1O M0. 

Such a direct connection of QSs with newly forming 
galaxies has been represented with models developing over 
the years from the first attemps by Cavaliere & Sza- 
lay (1986) and Efstathiou & Rees (1988), to Haehnelt & 
Rees (1993), Haiman & Loeb (1998), Haehnelt, Natarajan 
& Rees (1998), Cavaliere & Vittorini (1998, henceforth 
CV98). Concurring evidence is being provided by obser- 
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vations such as those by Fontana et al. (1998), and in 
Djorgowski (1998). 

Here we focus on the dramatic fall of the bright QS 
population later than z « 2.5, where more features can be 
observed and the picture is more complex. The time scale 
involved is around 2 Gyr, fast compared with the gentle 
demise of the galaxy formation envisaged by the hierar- 
chical clustering. The behavior of the optical luminosity 
function (LF) is dominated by luminosity evolution (LE, 
see Boyle, Shanks & Peterson 1988; Hatziminaoglou, van 
Waerheke & Mathez 1998), a trend also shared by the LFs 
in X-rays (Delia Ceca et al. 1994, Boyle et al. 1994) and 
in the radio band (Dunlop & Peacock 1990). 

Additional features appear in the optical data toward 
low z. First, a minor component of density evolution (DE) 
has been discerned at blue magnitudes Mb ~ —23 by 
Boyle et al. (1988), and confirmed by Maloney & Pet- 
rosian (1999). Second, La Franca & Cristiani (1997), and 
Goldschmidt & Miller (1998) found an excess of bright QSs 
over previous surveys; from these observations and those 
by Kohler et al. (1997) the LFs appear to be stretched 
out into a generally smoother shape on approaching low 
z. Finally, faint X-ray emissions and liner- type activity are 
increasingly found to be widespread among otherwise nor- 
mal galaxies, see Ho, Filippenko & Sargent (1997); Miyaji, 
Hasinger & Schmidt (1998); Giommi, Fiore & Pcrri (1998). 

On the other hand, sharp imaging of a number of bright 
QSs at intermediate z has revealed host galaxies which 
are not really newly forming, but rather are surrounded 
by a complex environment. Some still exhibit morpho- 
logical marks of a recent tidal interaction, in spite of the 
low surface brightness and short duration (~ 1/2 Gyr) of 
such features; others are found in close association with 
comparable or fainter companions; a fraction even appear 
to harbor secondary nuclei (see Hutchings & Neff 1992; 
Hutchings, Crampton & Johnson 1995; Rafanelli, Violato 
& Baruffolo 1995; Hasinger et al. 1997; Bahcall et al. 1997; 
Boyce et al. 1998; McLure et al. 1998; Ridgway et. al. 
1999). In addition, recent statistical evidence (Fisher et 
al. 1996; Jaeger, Fricke & Heidt 1998; Wold et al. 1999) 
points toward a poor group environment also for radio- 
quiet QSs, perhaps poorer than for the radio-loud objects 
(Yee & Green 1987, Lehnert et al. 1999). 

But the issue is not settled. Other observers detect 
at most marginal evidence that the environments of the 
bright QS are richer than the field (see Smith, Boyle & 
Maddox 1995; Teplitz, McLean & Malkan 1999; Croom 
& Shanks 1998), and find low incidence of peculiarities in 
galaxies active at the levels of the Seyfert nuclei (Malkan, 
Gorjian & Tarn 1998). Many point out the need to enlarge 
and complete the QS samples, and to compare them with 
proper controls. 

Such a complex data situation warrants exploring what 
predictions follow from pursuing the physical connection 
of bright QS activity with host interactions. Our guideline 
(taking up Cavaliere & Padovani 1989, Small & Blandford 
1992, CV98) will be that later than z w 2.5 the accre- 
tion feast ought to cease turning the QS rise into a de- 



scent, due to a gradual change in their environment. In 
fact, toward lower z major merging events (with the asso- 
ciated gas imports into the host) become rarer and rarer 
at galactic scales; the prevailing dynamical events are best 
described as interactions between developed galaxies, oc- 
curring mainly in the small groups which then begin to 
virialize. These events do not necessarily increase the dark 
halo masses, yet they can destabilize the gas in the hosts 
and trigger limited accretion episodes; these draw from 
the gas reservoirs in the hosts, and concur with starbursts 
in progressively exhausting them. But as exhaustion is 
approached other, lesser gas supplies become relevant for 
lower level activity, such as those provided from accretion 
of satellite galaxies that occours even in low density envi- 
ronments or in the field. 

We intend to single out and address directly the astro- 
physical processes involved, rather than base upon numer- 
ical and phenomenological treatments as in Kauffmann & 
Haehnelt (1999). Specifically, we shall simplify the dynam- 
ical history of the DM, to focus on the consumption and 
depletion of the gas reservoirs in the host galaxies. First, 
in §2 and §3 we will examine the outcomes of encounters 
of the host with sizeable companions, occurring mainly in 
poor groups. Then in §4 we examine accretion of satellite 
galaxies with their additional but small gas supplies, and 
consider also other processes that can occur in the field. 
Finally, in §5 we summarize and discuss our conclusions, 
stressing implications and predictions. 

2. QUASAR ACTIVITY TRIGGERED BY ENCOUNTERS IN A 

GROUP 

Our study will be based on the standard hierarchical de- 
velopment of cosmic structures past the era of galaxy for- 
mation; poor groups, in particular, are expected to start 
condensing at a considerable rate later than z « 2.5 in 
viable cosmologies. 1 

We hold that the approximate equality atzw 2.5 of the 
QS peak with the group onset constitutes not a coincidence 
but rather an intrisinc feature, telling of a connection of 
the QSs with the environment at scales larger than galac- 
tic. Our starting point is that QS hosts in early, poor but 
dense groups are bound to interact with their companions. 
The role of these dynamical events for nuclear activity is 
again to destabilize the gas in the host, specifically causing 
it to lose angular momentum. So they trigger inflows that 
rekindle the existing but dormant BHs into intermittent, 
dwindling episodes of supply-limited accretion and emis- 
sion. 

That encounters occur preferentially in groups is shown 
by many simulations (see Governato, Tozzi & Cavaliere 
1996, Athanassoula 1998). This may be anticipated (see 
CV98) on the basis of the average time r r ~ l/"g SI 1 ~ 
Gyr between encounters, basically because in groups the 
galaxy number density n g is high. On the other hand, in 
poor groups the velocity dispersion V only modestly ex- 
ceeds the internal galaxian velocities v g , so that the cross 
section for effective interactions (see Binney & Tremaine 
1987) is close to the geometrical value E ~ Tr(r g + r' g ) 2 

J For analytical evaluations we will use the critical cosmology with h = 0.5 and tilted cold DM initial perturbations; our numerical calcula- 
tions will also cover the flat, U = 0.3 case with h = 0.7 and standard cold DM perturbations (see Bunn & White 1997). We recall that in the 
former cosmology the standard hierarchical clustering implies the virializing masses to attain density contrasts close to 180 over the field and 
to scale up like M ~ 10 15 (1 + z)~ 6 /( n + 3 ) Mq, while the velocity dispersions scale like V oc M 1 / 3 (1 + z) 0,5 ; here n m — 2 is the power spectral 
index of the perturbations at group scales. 
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related to the radii r g of the host and r' g of the interaction 
partner. The time scale for inflows to develop is of order 
t w T g jvg ~ 10 _1 Gyr, the galaxian dynamical time. 

That such encounters can lead to substantial gas fun- 
nelling toward the nucleus is shown by a number of high 
resolution simulations comprising both the N-body treat- 
ment of the DM component and a hydrodynamic treat- 
ment of the gas (see Mihos & Hernquist 1996, Barnes & 
Hernquist 1998). These simulations show that over a few 
10 _1 Gyrs following a grazing encounter up to 50% of the 
gas orbiting in a galaxy at distances of kpc is driven into 
the central 10 2 pc, the resolution limit. The result fol- 
lows from a number of dynamical steps, as discussed by 
Mihos (1999). First comes the gravitational action from 
the partner, which depends on its mass M' and impact 
parameter b and lasts for a time b/V ~ r g /vQ w r. Then 
follows the response by the DM and the baryonic compo- 
nents developing in the host over times r; this is possibly 
amplified by internal metastability conditions to form true 
bars. Finally, it ensues the loss by the gas of a considerable 
amount of orbital energy and angular momentum \Aj/j\, 
which triggers the gas inflow. An internal loss of j (that 
is, Aj < 0) may be anticipated, since the asymmetries 
and inhomogencitics induced in the baryonic and in the 
DM components do not superpose or align; so the grav- 
itational coupling will transfer j from the former to the 
latter, more massive component. 

But most accretion episodes will be milder, originating 
from a fly-by of the host by a smaller or a comparable com- 
panion at impact parameters b ^ r g . From eq. (1.1) it is 
seen that in gas-rich hosts still containing m ~ 5 10 10 M Q 
of fractional gas mass Am/m <~ a few % sent to the 

central BH over At <~ t ~ 1CP 1 Gyr is enough to yield 
outputs approaching 5 1 46 ergs -1 . 

The actual gas amount made available by an encounter 
will limit the maximal luminosity attained by the QS be- 
fore fading out. The controlling factor will be |Aj/j|, 
which closely bounds from above the gas fraction / fun- 
neled to a small inward velocity (see Gunn 1977). In fact, 
the gas funneled inward may end up not only in accre- 
tion onto the central BH, but also in a less constrained 
nuclear starburst (see Sanders & Mirabel 1996), or it may 
be mostly dispersed. We shall assume that about 1/3 of 
the destabilized fraction / constitutes the fraction Am/m 
actually accreted. 

It is convenient to rewrite the output in eq. (1.1) as 



L <~ rj c 2 Ato/t oc m(t) f 



(2.1) 



apart from the factor 7/c 2 /3r fixed for a given host. This 
is to visualize and tell apart the long-term trend and the 
fluctuating component in the gas accretion. The factor 
m(t) drifts down on the long scale set by the time r r ~ 
Gyr between encounters; over the shorter scales r ~ 10 _1 
Gyr of a fly-by the factor m(t) stays nearly constant, but 
the other factor / 5; \Aj / j\ is statistically distributed fol- 
lowing the distributions of the dynamical parameters of 
the fly- by events. 

2.1. The shape of the LFs from encounter statistics 



Here we show how the distribution P(f) of the fractions 
/ gives rise to the shape of the LF, that we indicate with 
N(L). 

We expect that many encounters can destabilize gas 
fractions of order /{, w \Aj/j\ ~ 5%; these will lead 
to accretion of (Am/m)j ~ 2%, and to outputs up to 
Lb ~ 5 10 ergs -1 in a gas-rich host. We shall see that 
they produce LN(L) ~ const in this luminosity range. 
Higher luminosities require larger destabilized fractions 
/>/{,; in turn, these require closer encounters with larger 
companions, and these events are fewer. When we com- 
pute the corresponding distribution P(f), we expect it to 
drop steeply, producing a steep bright end of LN(L). 

To show this, we consider a group and inside it we exam- 
ine the interaction of a gas- rich host galaxy (with radius r g 
and circular velocity v g ) with a group companion of mass 
M' on its orbit with relative velocity V and impact param- 
eter b. The gas mass m in the host, at equilibrium on scales 
r of some kpes with the velocity v ~ (GMo/r) 1 / 2 enforced 
by the host mass M within r, has the specific angular 
momentum j — vr w GM a /v. The maximal Aj| obtains 
when the gas rotates slowly and responds strongly to the 
gravitational force from the companion. In such conditions 
(see Appendix A) we evaluate the time-integrated, relative 
change \Aj/j\ « j M' /M Q V b, which easily approaches 1. 
Using this to evaluate / in eq. (2.1), we see that L oc M'/b 
holds, on neglecting over the time scale r the slow decrease 
of m. 

So P(f) as defined above is computed on convolving 
the distributions of the dynamical parameters M' and b, 
namely, pi(M') and by P2(b) combined into 

P(f) = J JdM'db Pl (M')p 2 (b)S(f-AM'/b), (2.2) 

where the factor A = j /M D V is fixed for a given host and 
over times of order r. The distribution p\ (M 1 ) is dispersed 
over a wide range; we adopt the Press & Schechter (1974) 
form which includes the power-law factor M' _1 - 85 ) from 
M' = 510 10 M Q up to the mass of the host, often the 
brightest group member. 2 The scaling of P(f) in the 
upper range of / depends only moderately on P2(b). The 
latter reflects the poorly known geometry of the galaxian 
orbits in the group; but when the group is small, with 
membership of a few to several bright galaxies, the range 
of 6 extends by a factor of hardly 10, from the minimal 
value about r g (which sets encounters apart from prompt 
merging, see Governato el al. 1996) to the upper value 

RgN g (which defines the effective binary encounters 
in terms of the group radius Rq and the membership Afc)- 
We represent in fig. 1 the probability P(f) that we 
compute in the Appendix A from eq. (2.2) in three cases 
that bracket the variance expected. One is the simplest 
case A obtained on using a peaked distribution P2(b); we 
find that the shape in the main body of P(f) follows the 
shape of pi(M'), to yield P(f) = 0.85 (///b)" 1 - 85 ^" 1 be- 
fore going into a cutoff. The second case B is the other 
limit holding for a flat distribution P2(b), and yielding a 
shape still close to P(f) oc / _1 ' 85 before the cutoff; this 
will constitute our reference frame. Case C obtains when 
the gas rotation is fast and the gravitational torque from 



2 Similar results obtain on starting from a Schechter LF for the galaxies, see Zucca ct al. (1997), converted to mass distribution with the 
use of LocM 1±:L /3. 
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the partner can be described as a truly tidal effect (see 
Appendix A). Then integrating over the fly- by time we 
evaluate \Aj/j\ w jM'/M Q Vb x r/b; the shape of P(f) 
turns out to be somewhat flattened in the extended lower 
range of /. 




log f 

FIG. 1 — The probability distribution P(f) that a fraction / of 
gas is destabilized during an encounter, in the cases discussed in §2.1 
of the text: A (solid line); B (dashed line); C (dotted line). Press & 
Schcchtcr 1974 distribution for the partner masses from M' = 5 10 10 
to 510 12 M . 

To evaluate the corresponding shape of N(L), we con- 
sider that after times of order b/V ~ 10 _1 Gyr taken by a 
fly-by (see §2) the gas will increasingly flood the accretion 
disk in the nucleus on a galactic timescalc r again of or- 
der 10 _1 Gyr (see §2). The reactivated QSs will brighten 
up at the rate L = L/t, up to the value of L allowed by 
the accreted fraction Am/m //3; in the range L < Lb 
corresponding to minimal (Am/m)(, their number is con- 
served, that is, LN(L) ~ const applies. Then they fade 
off and drop out of the LF; the relative number of objects 
doing so in each range dL around L is given by the relative 
variation of P as given by 

d(NL)/NL = dP/P , (2.3) 

where the right hand side is evaluated at / oc L (note the 
analogy with the radiative transfer formalism) . The result 
is L N oc L^ 13 with /3 s=s 2 at luminosities brighter than Lb- 
In sum, past a break at Lb ~ 5 1 46 erg/s correspond- 
ing to / = ft, ~ 5%, the QS luminosity function fed by 
encounters steepens from N(L) oc L^ 1 to N(L) oc L^ 1 ^ 13 , 
and then goes into a cutoff. This evaluation will be checked 
with the systematic computation given in §3. 

2.2. Evolution from encounter rate 

The comoving LF evolves over the long time scales pro- 
vided by time between encounters r r = l/n g EV. The 



latter grows longer as z decreases, basically because the 
density n g oc (1 + z) 3 of galaxies in a group is propor- 
tional to the decreasing mass density in the background 
at virialization. This is offset only in part by the in- 
crease of the encounter velocities V oc Mq 3 (1 + z) 5 
caused by the gravitational clustering, which approxi- 
mately yields V oc (I+2) -3 / 2 in the critical universe where 
Mq oc (l + z)~ 6 holds. In the absence of strong galaxy evo- 
lution the cross section remains close to £ 47rr 2 «const 
in groups or poor clusters. Then t t oc (1 + z) -3 / 2 oc t 
closely obtains in groups. 

As a consequence, two components to the QS evolution 
may be anticipated. A first, milder component is due to 
the encounters becoming rarer as richer, but less dense 
groups come of age. This clearly causes in the QS popu- 
lation a mild DE proportional to t^ 1 oc (1 + z) 1 - 5 at late 
epoch z 1. 

The second and stronger component is in the form of 
LE; it occurs because on long time scales the luminosities 
decrease together with the available gas mass, following 
L oc m{t). In other words, the interactions also become 
less effective in feeding the BHs as the host gas reservoirs 
are depleted by all previous discharges. 

From P(f) given by eq. 2.2, we compute (/) = 
J df f P(f) « 15% corresponding to an average Am/m ~ 
5% and to a minimal fraction (Am/m)i ~ 2% actually 
accreted. So the gas is depleted following 

dm/dt w -(f)m/T r (t) , (2.4) 

to yield m(t) oc (t/t )-<f> W^", where t Q « 13 Gyr is 
the present epoch. The value r ro of the time between en- 
counters in local, virialized systems may be scaled from 
the classic census in the local field of bright galaxies with 
persisting signatures of interactions (Toomre 1977); this 
leads to a time around 2 10 2 Gyr between such encounters 
in the field. The value scales with the inverse of the galaxy 
density down to r ro « 1 Gyr in virialized groups, where 
density contrasts of about 2 10 2 over the field are indicated 
by the hierarchical clustering. The resulting behavior is a 
strong LE; in particular, in the critical universe the break 
luminosity Lb oc m(t) evolves as 

L b oc r 2 oc (1 + zf . (2.5) 

To sum up the contents of this Section, from direct eval- 
uation of the rates and strength of the accretion episodes 
driven by encounters we have found the following features: 
the comoving number of the QSs evolves like TV oc t^ 1 oc 
(1 + z) 1 - 5 for z < 1, the era of near completion for group 
formation, see §3; the shape of the LFs has the form of a 
double power law broken at Lb ~ 5 10 46 [(1 + z)/3.5] 3 ; the 
bright end is given by LN(L) oc (i/L fc )"' 3 with « 2. 
Collecting these results in one expression, we expect the 
LF of QSs to follow 

N{L, z) oc (1 + z) 1 - 5 L- 1 x {L/L b y p , (2.6) 

with the last factor applying only beyond the break Lb- 
These basic behaviors of N(L, z) will be checked and com- 
pleted in the next Section. 
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3. THE FULL LF FROM POPULATION KINETICS 

Here we present a systematic computation of N(L,t), 
the differential and evolving LF in comoving form, pro- 
duced when the BHs are reactivated by encounters to 
bright QS luminosities. We will base our computations 
upon the formalism of the so-called continuity equation 
(see Cavaliere et al. 1983). 

This pictures the statistics of the QSs as analogous to 
a hydrodynamical flow along the L-axis with sources and 
sinks. In fact, the equation that we are going to recall 
and use is the analog (whence the name) of the continuity 
equation for a 1-D fluid, with the density replaced by the 
comoving number of objects per unit luminosity N(L,t), 
and the velocity replaced by L. The equation stems from 
the general expression for the rate of change dt N in a given 
bin dL around L; the change is due to three causes. First, 
the objects are activated at the luminosity L and with the 
rate given by a source function S + (L,t). Second, they 
brighten up at the rate L on the time scale r discussed in 
§2, and so flow out of the bin dL at the net rate — 8l(LN). 
Third, they may quench for lack of fuel after attaining the 
luminosity L, at the rate given by a sink function S-(L, t). 

The algebraic sum of the above contributions yields the 
net rate of change dtN in the form of the kinetic equation 

d t N + d L (LN) = S+-S- . (3.1) 

This is the "continuity equation" treated by Cavaliere et 
al. 1983. We now discuss in turn the three astrophys- 
ical inputs L, S+, and S- that substantiate the general 
formalism for the present scope. 

In the present context, the brightening L = L/t de- 
velops on the time scale r ~ 1CT 1 Gyr (discussed in §2) 
taken by the gas to be directed on an inward course and 
to flood the accretion disk in the nucleus. This implies 
L = L i e( t ~ ti ^ T , where U is the epoch when reactivation 
and brightening begin from a faint luminosity L t . These 
"characteristic lines" L(t) are related to the streamlines 
in our hydrodynamical analogy, as shown in Appendix B; 
there we also show that they provide a simple way to find 
the relevant solution of eq. (3.1). 

The source function S+(L,t) represents the rate of BH 
reactivations per unit L. As to its L-dependence, reacti- 
vation implies scarce if any accretion occurring previously, 
so the initial luminosities Li are faint compared to the 
bright levels we focus on; hence to a first approximation 
we shall necglect the dispersion, and take S+(L) to be 
concentrated around a single value Li in the form of a 
5- function S(L — Li). As to the i-dependence, the reacti- 
vation rate is given by the (comoving) number density of 
hosts located in groups at the epoch t, divided by the mean 
time r r between reactivations. We focus on the member 
in a group that is so large and gas-rich as to fuel a sizeable 
activity upon encounters; often this will be single, based 
on the magnitude difference between the members of small 
groups observed by Ramella et al. (1999). So the number 
density of effective hosts will equal Na(t), the density of 
groups in the universe at the epoch t. In sum, our source 
function is given by 

S + (L,t)=N G (t)5(L-Li)/T r (t) . (3.2) 

Most rekindled sources attain L ~ Lb ~ 5 1 46 erg s _1 



as said in §2. So at faint luminosities L < Lb no quench- 
ing occurs, and only the source term S+ is present in eq. 
(3.1); then the solution of eq. (3.1) takes the simple form 
(derived in Appendix B, see eq. B.6) 

N{L,t) = T -^-L- 1 {L/L i )^ . (3.3) 

T r (t) 

Here we see the basic shape N(L) oc L^ 1 (corresponding 
to L N « const, as expected from our preliminary discus- 
sion of §2.1); this reflects the object population established 
by the brightening along the lines L — Li e i -?-~ ti )l T . The 
modulating exponent 7 (evaluated in Appendix B) takes 
values close to 7 = r/tc ~ 0.2 at epochs close to to (i.e, 
at z ~ 2.5), when the group number rises briskly. This is 
because at such epochs the source function 5+ oc Na(t) 
increases appreciably over the time scale r, so it tends to 
populate the faint bins close to Li faster than the objects 
can flow out to brighter bins, and N(L) steepens some- 
what. At z ^ 1 instead, Nc(t) is on its slow demise; so 
the exponent 7 vanishes or becomes slightly negative. At 
these late epochs the i-dependence of the LF is dominated 
by the decreasing factor r I T 1 (t) oc t~ x which yields a weak 
and late DE; this is maximized in the critical universe 
where it follows N oc t^ 1 oc (1 + z) 15 . But we recall 
from eq. (2.6) that for z 2.5 a stronger LE of the form 
Lb oc (1 + z) 3 is driven by the depletion of gas accreted or 
condensed into stars following each interaction. 

For L > Lb also the sink function S-(L, t) intervenes at 
the right hand side of eq. (3.1); this represents the number 
of objects per unit time that brighten beyond Lb, but then 
quench and drop out of the LF at some larger L due to 
fuel limitations. 

If these independently affect the QS population at a 
given z, then 5_ oc N is to hold. On the other hand, the 
drop out rate is related to the probability P(f) oc f~@ dis- 
cussed in §2.1; clearly, a steeper slope (3 will corresponds 
to faster drop out rates compared with the scale r of the 
brightening. These considerations lead to 5_ = N (3/t\ in 
other words, r//3 is the QS lifetime at high luminosities 
before fading out. Since (3 = \dlnP/dlnf\, and recalling 
that L — L/t holds, the relation also writes 

S-(L,t) = N(L,t)\L dlnP/L dlnf\ . (3.4) 

This is just the same as given by eq. (2.3) under the form 
d(N L) /dL = N L dP/dL P. 

The solution for L > Lb when both S + and S- apply is 
derived in Appendix B, eqs. (B.7) and (B.8), and yields 
the bright end behavior 

N(L, t) = (L/L^ (L/Lb)-^ 1 . (3.5) 

T r (t) L b 

The full shape of the LFs is computed in Appendix B 
using eq. (B.7). Specifically, if for 5_ we use the represen- 
tation given by eq. (B.9) with (3 =const, then the solution 
is analytic and reads 

t Nq (t) (Lj/LbT 1 

1 ' ] r T (t)L h (L/L b y+-f + (L/L b y+f+0 ' 

(3.6) 
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which gratifyingly coincides with the empirical fitting for- 
mula used, e.g., by Boyle et al. 1988, Pei 1995, La Franca 
& Cristiani 1997, Shanks at al. 2000. Note that this goes 
over smoothly from the limiting form (3.3) at the faint 
end to the other limiting form (3.5) at the bright end. 
Note also that the natural normalization of the LF given 
by eq. (3.6) turns out to be a few % shining QSs per 
bright galaxy; this arises with no parameter adjustement 
from the factor r/r r ~ 10 _1 , and from the fraction up to 
40% of bright galaxies residing in groups with membership 
M G > 3 (Ramella et al. 1999) which yields N G ~ 0.4/3 
relative to the bright galaxies. 




_ 2 , 1 1 1 

-10 12 




-2 ' ' ' ' ' ' 

-10 12 

log VL* 

FIG. 2 - Optical LF from interactions of the BH hosts in groups 
and in the field; the solid lines show the sum (see §3 of the text). Ef- 
ficiency r) = 10~ 1 , bolometric correction for the blue band Kg = 10, 
L, = 10 45 erg/s. a) Critical cosmology with h = 0.5, tilted CDM 
perturbation spectrum, b) Same for Qo = 0.3, Q\ = 0.7, h = 0.65, 
CDM perturbations. Data points from Boyle et al. 1988, and La 
Franca & Cristiani 1997. 



In fig. 2 we represent the results of numerical computa- 
tions from eq. (B.7) where we use the detailed numerical 
function P(f) (with / = i Lt /nine 2 ) given in fig. 1, case 
B. All computations are performed both in the critical 
universe with Q = 1 (fig. 2a), and in the low density, 
flat cosmology with n o = 0.3, fi A = 0.7 (fig. 2b). These 
lead to comparable results; actually, the latter case ac- 
cords marginally better with the data for QSs down to 
blue luminosities Lb ~ 10 44 erg s _1 . 

4. FUELING PROCESSES IN THE FIELD 

In fig. 2 we have added to the bright end of the LFs two 
small contributions which arise even when the QS hosts 
are located in the field, due to the following processes. 

(i) BH fueling from encounters of their host with a 
substantial partner; these occur mainly in the large-scale 
structures which contain most galaxies and have contrasts 
of a few units over the general average as shown by redshift 
surveys (Ramella, private communication). Such galaxy 
densities are still lower by about 10~ 2 compared with those 
in virialized groups (see footnote ^'), so the field encoun- 
ters take place at the correspondingly lower rate 10~ 2 r^ 1 . 
On the other hand, they involve more hosts by a factor 
about A/g/0.40 « 7; this arises (see the discussion at the 
end of §3) because the most effective groups have mem- 
bership JVq 3 and comprise about 40% of the bright 
galaxies. Thus the contribution from field encounters is 
around 7 10 -2 relative to encounters in groups. 

(ii) A lesser contribution (around 5 %, as discussed by 
CV98) is due to BHs flaring up to Eddington luminosities 
in the minority of galaxies still formed after z w 2.5 by the 
few major merging events still occurring. 

These two contributions are small but decrease slowly, 
because field galaxies undergo rare encounters and their 
gas is used mainly by the slow quiescent star formation 
(see Guidcrdoni et al. 1998). So on approaching low z 
such contributions are bound to emerge at the bright end, 
and to exceed there the rapidly decreasing contribution 
due to encounters in groups (singled out by the dotted 
lines in fig. 2). 

On the other hand, nuclear activities fainter than Lb ~ 
10 44 erg s _1 are less demanding in terms of accretion, and 
may feed on a number of different and superposing pro- 
cesses. To complete the above picture, we examine here 
the minor host-satellite interactions which become relevant 
when the gas reservoirs in the hosts approach exhaustion. 
From these events we anticipate a considerable addition 
to the faint end of the LFs; our expectations are outlined 
below, and are all checked and expanded in Appendix C. 

Since also the host galaxies located in the field can ac- 
crete their satellites, many more BHs can be activated 
by accretion related to such events. But lower activity 
levels L < rj-i 10 45 erg/s will be so allowed. This is be- 
cause of two concurring reasons: first, the satellites masses 
are small, in the range M s ~ 10 9 — 5 10 10 Mq comple- 
mentary to that discussed in §2.1; second, the usable gas 
supplies are provided mainly by the satellite cores, as we 
shall see below, and these contain meager gas amounts 
m < 10~ 3 M s . Very sub-Eddington L will be produced in 
this way, since the central BHs in the hosts are already 
grown up. Another expected feature is that the distri- 
bution of the gas masses to be accreted from the satellites 
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should be as steep as, or even steeper than the distribution 
of the faint galaxy masses. In conclusion, the associated 
LFs will be steep, and their normalization will be high 
due to the many field galaxies also activated in this mode, 
though at faint levels. 

On the other hand, the event rate will decrease sharply 
as satellites are accreted out of the initial retinue and their 
residual number per host galaxy is depleted. This will 
yield a strong DE (in fact, one which is L-dependent as 
we shall see), but will cause no LE since the satellite gas 
supplies will be accreted either all or nothing. 

Finally, as the accretion rate diminishes, the radiative 
efficiency itself may decrease below r\ ~ 10 _1 ; moreover, 
the optical emission is expected to be especially weak and 
often drowned into the galactic light, while the associated 
emissions in X-rays will be stronger and also more eas- 
ily pinpointed. These features arise because regimes of 
ADAF or ADIOS kind (see the discussion by Blandford & 
Begclman 1998 and references therein) are to set in at low 
accretion rates such as those fed only by satellites. While 
it is much debated which regime actually prevails, both 
will yield a low overall efficiency in using and/or convert- 
ing into radiation the gas that reaches the accretion disk. 
In both cases the hard X-rays should dominate over the 
soft X-rays, even more over the optical emissions, see Di 
Matteo ct al. (1999). So at these low accretion rates the 
soft X-ray and the optical bands comprise a small fraction 
of the small bolometric luminosity actually radiated. 

Our detailed computations are given in Appendix C. 
Here we stress that both the gas supply carried by a satel- 
lite, and the gravitational perturbation affecting the host, 
will be effective for accretion only when a satellite or its 
residual core sinks deep into the host body and delivers 
there the gas to be accreted. Such sinking may be triggered 
by encounters, but more often will begin as orbital decay 
caused by dynamical friction; the formal time scale is given 
by t s = 3 (M s /5 10 10 M Q ) Gyr, see Binney & Tremaine 
(1987). 

However, the effective time scale actually depends on the 
current satellite mass M s (t); it has been recently stressed 
how thoroughly the sinking satellites will be peeled off 
along their way by tidal disruption and stripping, espe- 
cially those which are initially larger, see Walker, Mihos & 
Hernquist (1996) and Colpi et al. (1999). The full process 
is very delicate to calculate and even to simulate numer- 
ically, but it basically leads to a weak dependence on the 
initial M s . In view of these uncertainties, our scope here 
will be limited to explore trends; to this specific purpose 
we adopt the phenomenological scaling 

r s = 3(M s /510 10 M Q )-« Gyr (4.1) 

with the exponent £ < 1 designed to describe a weak de- 
pendence on the initial M s . 

For consistency, we have also to consider that the gas 
effectively sunk into the host body is provided only by the 
fraction residing in the satellite cores; this is not only small 
but also will depend weakly on the initial total mass M s . 
Again phenomenologically, we will use the exponent 8 < 1 
to describe this dependence in the form 

m = 510 7 M Q (M S /510 10 M Q ) 9 Gyr . (4.2) 

The full LFs so generated are computed in Appendix 
C, and given in eq. (C.5), using the appropriate form of 



the kinetic equation (3.1). To illustrate the generic behav- 
iors we give the result corresponding to the specific values 
£ = 6 = 0.7, which yield 

N(L,t) cx L- 2 - 2 e - Lt/3L ° . (4.3) 

Here L s ~ 10 45 erg s~ x is the maximal bolometric lumi- 
nosity from the gas mass m = 5 1O 7 M accreted from the 
core of a large satellite which spirals into the host on the 
scale t s w 3 Gyr; the time t is also given in Gyrs. 

6 I 1 




- 
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log L x /L t 

FIG. 1 - The contribution (dotted lines) to the X-ray LFs in the 
keV range by accretion of satellites with masses below 510 10 M Q ; 
we have taken i|/kj = 10~ 3 as discussed in §4 of the text. This 
contribution adds as shown by the solid lines to the one due to en- 
counters in groups, represented in fig. 2a. Data points from Miyaji 
ct al. 1998. 

In fig. 3 we represent specifically the soft X-ray LFs at 
low z resulting in the critical universe from cannibalism of 
satellites. We have adopted tj/kx ~ 10~ 3 , considering (as 
commented above) the small efficiency r) < 1CP 1 and the 
large bolometric correction nx > 10 expected for soft X- 
ray emission at low accretion rates. We represent also how 
these LFs add to the contribution from group encounters 
that has been computed in §3 and represented in fig. 2a. 

In sum, the contributions from satellite cannibalism to 
the LFs tend to pile up toward the faint end, producing 
steep LFs which steepen yet to a cutoff. Such contribu- 
tions will evolve following a strong, L-dependent DE. 

5. CONCLUSIONS AND DISCUSSION 

We have shown that a steep QS fall for z < 3 is to re- 
sult from the effects of the hierarchical growth of structures 
around accreting BHs, amplified by the exhaustion of the 
gas available for the accretion. 

The hierarchical clustering envisages the era of galaxy 
formation to go over to an era of group assemblage for 
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z 2.5. In a group, the settled host galaxies keep on 
evolving for a while, no longer by major merging events but 
rather via encounters which occur at rates t~ x ~ n g E V ~ 
Gyr" 1 . 

These are milder events, yet able to destabilize the host 
gas, specifically by causing considerable loss of angular 
momentum as tackled in §2. We have shown how they 
trigger intermittent and supply-limited accretion episodes 
onto the central BH, yielding blue luminosities Lb ~ 10 45 
erg/s for some 10 _1 Gyr in large, gas-rich galaxies. But we 
have additionally shown in §2.2 how such substantial gas 
discharges, which also fuel simultaneous starbursts activ- 
ity, deplete the host reservoirs and cause the gaseous mass 
to decrease substantially following m(t) oc t~ 2 on average. 

Thus the encounters not only turn rarer (r r oc t) as the 
clustering proceeds toward groups richer but less dense, 
and so produce a later, weak DE; but they also become 
less and less effective in terms of triggered accretion, and 
this causes strong LE with the luminosities scaling down 
close to Lfc oc (1 + z) 3 . The result we predict for bright 
QSs in the range z ^ 2.5 is outlined in eq. (2.6), is derived 
and shown in more details in §3 to yield eq. (3.6), and is 
represented in full by fig. 2. This agrees not only with 
the observations by Boyle et al. 1988, but also with the 
recent re-analysis of variuos samples by Maloney & Pet- 
rosian 1999, and with the greatly extended samples being 
provided at M B < -23 by Shanks at al. 2000. 

As to the shape of the LFs, we have shown in §2.1 that 
the statistics of the encounters naturally yields a double 
power-law broken from about L~ 12 toward the faint end, 
to L~ 3,2 and steeper at the bright end, in accord with the 
above data. But toward low z a smoother overall shape 
is produced by the changing tilt 7 given by eqs. (3.6); 
this concurs with the enhancement of the bright end at 
Lb ~ 10 46 erg/s provided by the additional feeding pro- 
cesses in the field discussed in §4 and shown in fig. 2. Our 
result accords with the observations by La Franca & Cris- 
tiani (1997) and by Goldschmidt & Miller (1998), which 
are confirmed by a recent survey being analyzed by Cris- 
tiani and collaborators (see Grazian et al. 2000). 

The above values and behaviors are given for the criti- 
cal universe, see fig. 2a; but fig. 2b shows that similar or 
better results obtain in a flat, low density cosmology. 

From a different path, we agree with Cattaneo, Haehnelt 
& Rees 1999 and Kauffmann & Haehnelt 1999 that strong 
QS evolution down to very low z implies destabilization 
of substantial fractions of the gas stockpiled in the hosts. 
Only a minor part of such gas ends up into accretion 
episodes; the major part is likely to be heated and/or 
ejected, or used up to produce - over and above the quies- 
cent star formation - star bursts specifically in the bright 
host galaxies, mostly ellipticals. But we see reasons why 
the LE evolution component may be expected to slow 
down on approaching z 0.5. For example, when interac- 
tions cause outright galaxy aggregations the host mass M g 
grows appreciably while its internal density p g decreases; 
then it is easily checked that the scaling of t t includes the 
factors Mg^ 3 p 2 J z giving rise to t t oc t 1+t with e < 1. Cor- 
respondingly, from eq. (2.4) the late evolution L oc m(z) 
brakes into a softening exponential at low L and z. 

On the other hand, when it comes to faint AGNs at 
z ^ 0.5 the situation becomes more complex, since the 



underlying BHs may be fed at low levels by various pro- 
cesses which superpose. We have explored in §4 and in fig. 
3 one such process, namely, the accretion of satellites with 
the associated gas supplies; this contributes many more 
but weaker sources, in the range of the Seyfert galaxies 
and below. These sources undergo a strong L-dcpcndcnt 
DE, similar to that observed in X-rays by Miyaji, Hasinger 
& Schmidt (1998); many, in fact, will be more easily pin- 
pointed in X-rays at Lx ~ 10 43 erg/s or fainter. 

Other feeding processes are conceivably provided by the 
self-destruction of giant nuclear star clusters (see Norman 
& Scoville 1988), and by internal instabilities and bar for- 
mation in disk galaxies (see Heller & Shlosman 1994, Mer- 
ritt 1998, Scllwood & Moore 1998). It will be interesting 
to see what specific shape and evolution are predicted for 
the LFs of AGNs fed in this latter fashion. 

Here we focussed on dynamic losses of the gas angular 
momentum j as primary triggers of the inflow eventually 
leading to BH accretion; the gas may undergo other ad- 
ventures before actually fueling the central BH, see e.g. 
Sicmiginowska & Elvis (1997). But the straight prediction 
from our study substantiating the outline by Cavaliere & 
Padovani (1989), is that the emissions from the average 
QS or AGN should go progressively sub-Eddington later 
than z w 2.5, even at constant r\ <~ 10~ 4 ; this accords 
with the data trend taking shape from various observa- 
tional lines, see Sun & Malkan (1989), Haiman & Mcnou 
(1998), Salucci et al. (1998), Wandel (1998). 

Our specific prediction for the average Eddington ratio 

is 

L/L E » [(1 + z)/3.5] 3 M BH (2.5)/Mbh(z) , (5.1) 

which goes down to about 10~ 2 by z 0.1. This is based 
on the intrinsically depleted accretion onto already exist- 
ing, massive BHs which yields L oc (1 + z) 3 as discussed 
above; the other factor is the related mass growth by the 
ratio Mbh(z)/M B h(2.5), see CV98. 

In fact, as outlined there and expanded elsewhere, by 
all episodes of rekindled accretion in the range z ^ 2.5 
the central BHs now dormant in most galaxies grow by 
an average factor Mbh(0)/Mbh(2-5) ~ 5. This will be 
lower for hosts which are disk dominated at radii of a few 
kpes, based on the values of Aj/j indicated in §2.1 by 
case C compared with case B. Otherwise, the factor will 
be limited by the bound Mbh ~ 10~ 2 Mbuige set when 
a large central condensation enforces axial symmetry de- 
spite the dynamical perturbations. Correspondingly, the 
local density of baryons collapsed into BHs amounts to 
Pbh ~ 310 15 /i 3 M Q Gpc -3 . These findings agree with 
the observational estimates of the central BH masses in 
many local galaxies given by Richstonc et al. 1998, and 
with the upper bounds set by Magorrian et al. 1998. 

Toward a comprehensive picture, we point out that the 
losses of j during galaxy-galaxy interactions as here ana- 
lyzed go back to dynamically generated asymmetries of the 
gravitational potential. Compared with the major merg- 
ing events during the build up of the spheroids, these do 
not differ in kind, but only as for the lesser strength and as 
for their falling rate and decreasing gas amounts involved. 

So we conclude that most of the bright QS evolution 
can be understood in terms of one engine (the accreting 
BH) and one trigger (dynamical destabilization of the gas), 
which however operate under two adjoining regimes geared 
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to the development of cosmic structures. Such regimes are 
related to the early, violent merging events when the gas- 
rich spheroids form; and to the later, dwindling interac- 
tions of host galaxies in groups or even in the field, with 
the slow demise of galaxy formation magnified by the fast 
depletion of the gas supplies. 

We thank for helpful exchanges B. Guiderdoni, B. Lib- 
erty P. Rafanelli and M. Ramella. Early stimulating dis- 
cussions with G. Hasinger and C. Norman are gratefully 
ackowledged. Thanks are due to our referee for stimu- 
lating us to clarify our presentation throughout the MS. 
Work supported by partial grants from ASI and MURST. 

APPENDIX A 

As said in §2.1 of the main text, we consider in a group 
of radius R G and membership Af G > 3 a gas-rich host 
galaxy with radius r g and circular velocity v g , and a com- 
panion of mass M' flying by with impact parameter b and 
relative velocity V. The gas in the host has specific an- 
gular momentum j = vr « GM Q /v if it is at equilibrium 
on scales r ~ kpc with the velocity v w {GM /r) x l 2 under 
the gravity provided by the host mass M Q within r. 

Let us consider a ring of gas. When the rotation is slow 
with v/V < r/b the gravitational torque exerted by M' on 
the unit mass may be evaluated from the angular deriva- 
tive of the interaction energy 



factor r/b, so that 



dE v 2GM'r 
dd> irb 2 



(Al) 



The overall variation Aj of the specific angular momen- 
tum will result from the torque time-integrated along the 
partner orbit: 



dtT 



GM'r 



Vb 



to read 



M'j 

MnVb 



(A2) 



(A3) 



Approximately 1/2 of the gas is braked and flows toward 
the nucleus, while the rest flows outward. 

With accretion so triggered, the probability P(f) of 
destabilizing a fraction / = AM' jb before fading out is 
given by eq. (2.2) in the main text. In the limiting case 
(denoted by A in the text) of a peaked distribution P2(b), 
the integrations are easy to perform analytically, to find 
P(f) = 0.85 {f/h)~ fb 1 ' havm g assumed the Press & 
Schechter shape pi(M') oc M' -1 - 85 . This result is repre- 
sented in fig. 1 by the solid line. 

The opposite limit, case B, obtains on considering a flat 
distribution P2(b) in the range r g - R G M G 1/3 . Then 
the numerical integration of eq. (2.2) yields the dashed 
line in fig. 1. This accords with the analytic evaluation 
which provides the shape P(f) oc (///b) -1 ' 85 between 

the limits AM' mm /r g < / < AM' max R G Af G 1/3 , before 
the final cutoff. We have used M' min = 5 10 10 M and 
M' max = 5 10 12 Mq, as said in §2.1 of the main text. 

Case C is one of relatively fast rotation with v/V > r/b. 
Then the net torque may be computed from differentiat- 
ing eq. (A.l) across the diameter. This brings in the tidal 



M'jr 
M Vb 2 



(A4) 



now holds. The luminosities so produced are given by 
L = r)c 2 m(t) M' jr/M Vb 2 r. Note that jr/M r is 
nearly independent of the host parameters; this is seen 
on recalling that v 2 « GM /r, and using the approxi- 
mation r « r/v. The main dependences of L are on m, 
which however stays nearly constant over the time scale 
of a fly-by, and on the interaction parameters which are 
statistically dispersed as here computed. 

For the purpose of the convolution (2.2) we now have 
/ = Ar M'/b 2 , and as before we have to integrate over 
the distributions of M' and b. Considering again a flat 
distribution for P2(b) we numerically obtain the dotted 
line in fig. 1. The analytic evaluation yields again a 
shape close to P(f) oc / _1 ' 85 , but now in the range 

ArM' ml Jr 2 g < f < Ar M' ma J 'R 2 G M G 2/S this is followed 
by a cutoff, and is preceded by a somewhat flatter section 
toward the extended lower end. 

APPENDIX B 

Here we give the analytic solution of the kinetic equation 
cq. (3.1), namely of 



d t N + d L (LN) =S+-S-, 



(B.l) 



where specifically L = L/t applies as explained in §3 of 
the main text. 

The simplest presentation is based on rewriting the 
2 nd term on the left hand side in the form 8l{LN) = 
L 8lN + N 8lL; so we may look at the equivalent equa- 
tion 



d t N + L d L N = -N d L L + S+-S- 



(B.2) 



The advantage is to visualize that the two terms on the 
left hand side make up the differential operator dt + L 
this is the total time derivative along the lines (known as 
"characteristics", see Shu 1992) defined by the solutions 
L = Lie( t ~ ti " , / T of the equation L = L/t, where U is the 
epoch when reactivation and brightening begin from the 
luminosity Lj. Such a derivative, that we denote with a 
dot, corresponds just to the Lagrangian derivative along 
the streamlines in the hydrodynamical analogy presented 
in §3. 

Thus along the characteristics we can reduce our partial 
differential equation eq. (B.l) to the ordinary differential 
equation 

N = -N^r + S+ - S- . (B.S) 
dL 

This is l s *-order, linear and in canonical form, with solu- 
tions given in textbooks, see Dwight 1961. 

We are interested in the initial condition N(L,t) = for 
t — > 0; that is, very few QSs were activated by encounters 
in groups at early epochs corresponding to z > 2.5 when 
few groups existed. We first write the solution in the faint 
range L < Lb where 5_ = holds, as explained in §3; 
there we have 

N(L,t) = [ dt' S+(L',t')e-ft> dt " dL/dL . (BA) 
Jo 



10 



Actually this expression simplifies a great deal because of 
two circumstances: first, dt = dL/L holds on a stream- 
line, so that we may equivalently integrate over dL; sec- 
ond, the expression of S + eq. (3.2) is proportional to 
S(L — L^. These two features together make the inte- 
grations straightforward, and the result is 



N(L,t) = 



N G (U) 

T r (ti)L 



(5.5) 



To make explicit the contents of this solution we sub- 
stitute L = L/t, and express ij in terms of t using 
t — ti = rlnL/Li along the characteristics. The differ- 
ence t — ti matters close to the epoch of group forma- 
tion t G (that is, z w 2.5) when the group number density 
rises rapidly, following closely N G (t) — N G (ti) e(*~* l )/* G ; 
so N G (ti) = N G (t){L/Li)- T l tG may be substituted in eq. 
(B.5), while T r (tj) « T r (t) holds. Thus the explicit solution 
reads 

TN G (t) 



N(L,t) = 



T r {t) 



'-L- 1 (L/LiY 



(5.6) 



where 7 = r/t Gl This is the eq. (3.3) of the main text. 
We end the study of the range L < Lb by noting that the 
integrand —dL/dL in the exponential of eq. (B.4) is just 
the coefficient of the term proportional to N in eq. (B.3). 

When we consider the bright range L > Lb, we 
see from eq. (3.4) that also the sink function 5_ = 
— N L dlnP/L dlnL intervenes. Now the complete co- 
efficient of N reads —dL/dL + L dlnP/L dlnL, and this 
will appear in the exponential of the solution which now 
writes 



APPENDIX C 

Here we compute in detail the contribution to the LFs 
from accretion of satellite galaxies, following up the dis- 
cussion in §4. 

We consider satellites of dynamical masses in the range 
M s - 10 9 - 5 10 10 Mq, with the Press & Schechtcr (1974) 
distribution N S (M S ) oc M" 1 - 85 . We assume a satellite to 
contain a usable gas mass m ~ 10 _3 M S , since only the 
gas residing in the core will be ultimately effective for BH 
feeding. 

We use the same basic formalism of the continuity equa- 
tion for N{L,t) as given in §3 of the main text by eq. (3.1), 
with L — L /r. But here the form of source S+ and of the 
sink function S- will require a number of changes. They 
go back to the following facts: here all dormant BHs in 
bright galaxies, whether in groups or in the field, may be 
activated; in addition, Am ~ m holds (the gas in the 
satellite cores is accreted all or nothing); finally, accreting 
these meager gas supplies yields relatively low bolometric 
outputs L ~ r/mc 2 /r bounded by the value L s ~ 10 45 
erg s _1 for a large satellite. 

The source function S + is now proportional to the 
number density Nbh of all BHs formed by z = 2.5 
which amount to about 1 per bright galaxy (see Hachnclt, 
Natarajan & Rees 1998); but it is also proportional to the 
diminishing satellite number per host galaxy J\f s if) • These 
factors are divided by the time r s taken to completely ac- 
crete a satellite, see eq. (4.1). The result reads 



S + = N BH N s {t)5{L-L i )/T s 



(CI) 



N(L,t) = / dt' S+{L',t')e 
Jo 

(5.7) 

A similar calculation applies as that described for eq. 
(B.4); the term added to the exponential will give rise 
to a factor multiplying the faint end form given in eq. (B. 
5). Using the expression P(L) oc (L/Lb)^ 13 (0 « 2) found 
in 52, the result reads 



•/*<**" (dL/dL-LdlnP/L dlnL) In turn, the satellite number Af s (t) is governed by 

dMs/dt = —J\f s (t)/r s , assuming the satellites are indepen- 
dently accreted. This yields the long term decrease 

N s {t) = Ki e-t'T' (C.2) 



N(L,t) = 



tN g (U) 
T r (ti) L Lb 



(B.S) 



from the initial retinue N S i that we conservatively take at 

Ki - 10. 

The form of the sink function S- again follows eq. (3.4), 
namely 



Recalling that tj may be expressed in terms of t using 
ti = t — rlnL/Li on the characteristics, we obtain eq. 
(3.5) of the main text. 

To compute the full LFs, we use for the coefficient 
S-(L,t)/N (which in fact constitutes the integrand in the 
second exponential of eq. B.7) a smooth transition across 
the break Lb- This start from zero at L <C Lb where S*_ 
vanishes (see §3 and above), and for L > Lb has to go to 
the full value given by eq. (3.4). A convenient representa- 
tion is provided by 



S-(L,t) = N(L,t)\L dlnP/dL\ . 



(C.3) 



S_/N = (3{L/Lbf/T[l + (L/L b f] 



(5.9) 



But now P(L) dL = N{m) dm is given straighforwardly in 
terms of the distribution of the satellite gas masses, since 
L oc Am « m hold. As explained in §4 of the main text, 
we have to consider that the effective gas supply is pro- 
vided by the satellite cores. Then the effective gas mass 
will depend only weakly on the total satellite mass, fol- 
lowing m oc Mf (with 9 < 1) as indicated by eq. (4.2) 
and related comments. So N(m) dm oc M" 1 - 85 dM s holds, 
and P(L) oc (L/L b )~ 1 ~ - 85 / 6 ' obtains for bolometric lumi- 



nosities exceeding Lb 



10 4 



corresponding to the 



This may be used to integrate eq. (B.7); if we keep 
=const, we obtain the analytic solution given by eq. 
(3.6), which has the form of a widely used fitting formula 
for the LFs (see text at the end of §3) . 

Alternatively, we may integrate numerically eq. (B.7) 
using the numerical form of P(f) given in fig. 1, to com- 
pute the LFs represented in fig. 2 of the text. 



lower end of the satellite dynamical masses. 

This we use in the solution of the kinetic eq. (3.1), al- 
ready given in general form by eq. (B.7). Evaluating the 
integrals as described there, we find 



N(L,t) 



T N B flNs(t) (£/£ fe )-2-O.85/0 



(CA) 
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But following eq. (4.1) of the main text and related com- 
ments, we have to consider that t s oc (with £ < 1) 
holds together with the dependence lama used 
above; it follows that t s oc L~^/ e . When we make ex- 
plicit this dependence and the decrease of Af(t) given by 
eq. (C.2) the result reads 

N(L,t) = TN ™ Msi (L b /L s y+°- S5 /<> x 

(L/L s )- 2 +^-°-^/ e e -(Vi s ) s/e */3 . (C . 5) 



Here all times are given in Gyrs; L s ~ 10 45 erg s" 1 is the 
maximal bolometric luminosity activated by the accretion 
of gas mass m = 5 10 7 M© from the core of a large satellite 
of mass M s = 5 10 10 M & spiralling down on the time scale 
t s = 3 Gyr (see the coefficients in the eqs. 4.1 and 4.2). 

Note that the QSs or AGNs fueled by satellites outnun- 
ber those activated by encounters in groups by the factor 
NBt{N{t) T r /NcT s ~ 5; this specifies the high normaliza- 
tion anticipated in §4. Eq. (4.3) in the main text obtains 
for the illustrative values £ = 8 = 0.7. 
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